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ABSTRACT 

£\J ' We briefly review the current situation with various relations between knot/braid polynomials (Chern-Simons correlation 

functions), ordinary and extended, considered as functions of the representation and of the knot topology. These include linear 
skein relations, quadratic Plucker relations, as well as "differential" and (quantum) *4-polynomial structures. We pay a special 
attention to identity between the yi-polynomial equations for knots and Baxter equations for quantum relativistic integrable systems, 
related through Seiberg-Witten theory to 5d super- Yang-Mills models and through the AGT relation to the g-Virasoro algebra. 
This identity is an important ingredient of emerging a 3d — bd generalization of the AGT relation. The shape of the Baxter equation 
(including the values of coefficients) depend on the choice of the knot/braid. Thus, like the case of KP integrability, where (some, so 
far torus) knots parameterize particular points of the Universal Grassmannian, in this relation they parameterize particular points 
in the moduli space of many-body integrable systems of relativistic type. 



p ■ 1 Ordinary and extended HOMFLY polynomials 

The HOMFLY polynomials [1] 

H%( q \A)\ A=gN = (Tr R Pe, P £A) su{N) (1) 



X 



are the analytical continuation in N of the Wilson-loop averages in 3d Chern-Simons theory [2] with the gauge 



group SU (N) and the action S = j d3x Tr (AdA+^A 3 ) , where q = exp-^jj. They can be lifted to extended 

polynomials [6, 7]: %j|{<z|Pfc}, which depend on infinitely many time-variables {pk} and on particular braid 
representation B of the oriented knot K (the braid arises in projection onto a 2d plane and depends on the 
1 particular embedding of the knot and on the choice of the plane), so that 



H*(q\A) = {q-^ R A-\ R \f {B) nl{q\ P t} (2) 

where is the writhe number, i.e. the algebraic number of crossings in the braid B, and 

A k _ A -k 

Pi = — — 3 

q K — q K 

At this topological locus the r.h.s. of (2) does not depend on the choice of B for a given /C. 

An m-strand braid is parameterized by a sequence of integers with i = l,2,...,m — 1 and j — 1, . . . , n, 
their meaning can be understood from the picture (in this figure a\ = —2, b± = 2, 0,2 = — 1, &2 = 3: this is knot 
810): 




Representation is parameterized by the Young diagram: an ordered set of l(R) positive integers R — {r% > 
Ti > ■••> r l{R) > 0} = [n,r 2 ,...,n (fl )]. 

One of the hot questions in Chern-Simons theory is to find and investigate the equations, which the HOMFLY 
and extended HOMFLY polynomials satisfy as functions of cijj and r^. 
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The main source of information about these relations (besides pure empiricism) is the character expansion 
separating dependencies on the braid and the time variables 1 : 

n%{q\ Pk }= ]T c B RQ (q)S Q { Pk } (4) 

Qhm\R\ 

where Sq are the Schur functions, i.e. the characters of GL(oo) in representations (Young diagrams) Q of the 
size |Q| = m\R\ and c R g(q) are represented as traces in the spaces of intertwining operators of combinations of 
simple matrices, 

4 Q (q)^Tr Q ]J^J[n^u}j (5) 

described in detail in refs.[6, 7, 8, 9]. More artful character expansions (especially in the Hall-Littlewood or 
MacDonald polynomials) can be even more useful for particular applications [10, 11, 12]. 

In this paper we briefly review the main known linear and bilinear relations between the HOMFLY poly- 
nomials, as well as some closely related topics. The subject is still at empirical stage of development, far 
less understood than, say, Virasoro constraints, various recursions and integrability in matrix models, though 
a list of known properties and facts is quite impressive. Therefore, it deserve an attention and further deep 
investigation. 



2 Skein relations for ordinary and colored HOMFLY polynomials 

The oldest known equation of this type is the skein relation, which for extended polynomials states that 

%(••• •••) _ %(••• o«-l •••) = [q- q- 1 ) % { n - a ^ - ) (6) 

for any particular parameter ay. The origin of this relation is that the 7£-matrix, acting in the channel 
[1] <E) [1] = [2] © [11] has two eigenvalues, q and —l/q, associated with the two irreducible representations [2] and 
JTl] respectively, and thus satisfies the Hecke algebra constraint 

(ll-q^K + q- 1 ) =0 (7) 

Eq.(6) is then a direct corollary of representation (5). 

For bigger representations R, the product i? <g> R — (BsS and 7?.-matrix in this channel has many different 
eigenvalues: as many as there are irreducible representations S in this expansion. These eigenvalues are equal 
to esq* s ', where ks is an eigenvalue of the cut-and-join operator [13]: W[ 2 ], associated with its eigenfunction 
character Ss{pk}, an d es = ±1, depending on R and S. Therefore, this 7\L-matrix satisfies 

11 (ll- esq" 3 ) =0 (8) 



and, as a corollary of (5), 



for any variable ay . 

For example, in the case of the first symmetric representation, one has 

(n - q 6 ) (n + <? 2 ) (n - 1) = o (10) 

i.e. 

U[ 2] aii+3 - - (1 + q 2 + q 6 )H^ - ) + q 2 (l + q 4 + q 6 )H^ ^ - q 8 H[ 2] ^ ^ = (11) 

while for the antisymmetric representation 

a w+ 3 ...) _ (1 + g -2 + q-^u^ ^ -) + g -2 (1 + q -A + q-^H^ ^ +1 ^ - q^H^ °" = (12) 



1 This presentation is essentially based on earlier attempts in [3, 4, 5]. 
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To make use of the skein relations for evaluation of the HOMFLY polynomials one needs to supplement them 
by "the initial conditions": in the case of R = □ one can express arbitrary V. aij with arbitrary {a^} through 
those where all = or 1. These polynomials should be found by some other means (e.g., at topological 
locus topological invariance can be used to decrease the number of strands). For higher representations the 
"boundary" ay can take more values. 

For superpolynomials [14, 15] the skein relations are non-trivially deformed. Already for the simplest 2-strand 
torus knots and links only one half of the relations remain intact: 

p[2,n+l] _ p[2,n-l] = (q - q-^-pi 2 '™] (13) 

for odd n (so that there are links at one of the sides of the relation), but for even n this is not true. The origin 
of this phenomenon is that the character expansion (in MacDonald polynomials) is different for knots and links 
[10]: 

f - (1+ fJ q % e) ■t- n M [11] {p} for n odd (knots) 
7?P."1 = q n M( 2] {p} + } (14) 

{ + g (i-" g 0) -t- n M {11] {p} for n even (links) 

3 Other relations for knot-dependence 

Skein relations restrict the dependence of the knot polynomial on the shape of the braid. They are the only 
universal equations of this kind, known so far. However, it is already clear that they are not the only ones: 
some traces of a more profound structure are already observed. We give just two examples, both for torus knots 
and links. 

3.1 "Evolution" of knots induced by skein relations 
Torus knots/links 

Skein relation as it is can not be directly deals with only torus knots, except for the 2-strand case: for the 
number of strands to > 2 it mixes torus knots with non-torus ones. However, the simple procedure, described 
in the previous section, immediately provides another kind of relations which provides a set of skein relations 
closed on only the m-strand torus knots. Namely, the HOMFLY polynomial for the m-strand torus knot/link 
[to, n] is a weighted trace of the n-the power of a special combination of m — 1 "/^-matrices, 

?4"'" ] M = Tr7C, K m =K 12 K 2Z ...K m ^ m (15) 

If instead of an individual 1Z in the standard skein relation, one now diagonalizes TZ m , one can immediately 
write the analogues of eqs.(8) and (9) for the series of knots [m, n] with given m. For example, for the 3-strand 
torus links in the fundamental representation one has [7] 

U [ n' n] {p] = q 2n S [3] {p} + Tr 2x2 (KUKWy •5 [21] M + (-l/ <Z ) 2 "5 [111] M (16) 

and the 2x2 matrix WJVJJ^ has eigenvalues e ±27 ™/ 3 . This means that our 1Z 3 in this case has four different 
eigenvalues: q 2 , e ±27Tl / 3 , q~ 2 , i.e. satisfies 

(n 3 - q 2 ) (n 3 - q - 2 ) (n 2 + n 3 + 1) = o (17) 

and this implies the equation for torus HOMFLY polynomial: 

(ed/On _ ^ ^a/On _ q -2^ ^28/dn + p a/8n + {p} = 

= {nt l+4] {p} + n^ n+3] {p} + nt l+2] { P }) 
-d 2 + T 2 ) (nt l+3] { P } + n l ^ n+2] {p} + nt l+1] { P }) + 

+ (n& n+2] {p} + H^ n+1] {p} + H& n] {p})=Q (18) 
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Thus one obtains a kind of evolution in the number of TZ m . This evolution can be described by "an evolution 
operator" [10], which is nothing but the cut-and join operator [13]. Moreover, this evolution can be naturally 
continued to the superpolynomials [10, 11]. 

If one wants instead an equation for the topological invariant HOMFLY polynomial, it is necessary to restrict 
time- variables to topological locus (3), pu — p* k and also include the proper normalization factor (2). In this 
case the evolution can be considered as an evolution in the number of strands. 



Non-torus knots/links 

In fact, in just the same way one can write a version of skein relation for an arbitrary one-parametric family 
of braids, obtained by iteration of any braid group element: B^ = BoB™. If the product 1Z\ of 7£-matrices, 
associated with B\, has eigenvalues then 

nfa-&(«))=° 

i 

and the corresponding extended HOMFLY polynomials satisfy 

U{ ed/dn -ZM)K BM {p} = (20) 

i 

independently of the choice of the "initial" braid element Bq. Thus, one can start an evolution with any "initial 
condition" Bo. For examples of such sequences see [10, 7], note that eigenvalues can be somewhat non-trivial 
functions of q. 

Generalization of this procedure to superpolynomials seems straightforward, but it is still unclear if it really 
works in full generality, see [10]. 



3.2 Strand-changing relations 

Like the skein relation itself, the modification, considered in the previous subsection does not change the number 
of strands in the braid. However, the strand-changing is also possible. An example is the following relation, which 
at the moment is very restricted: only torus knots, only SU(2) gauge group (A = q 2 , i.e. Jones polynomials), 
only fundamental representation, but instead it is in a very different direction in the braid-parameter space, 
and it survives (with minor modifications) a /3-dcformation to superpolynomials [16]. 

Namely, for torus knots of the special type [m, m+ 1] (m-strand braids, but a specially restricted n = m+ 1) 

n^- m+1] {p} = (T'-^hM - <? (ro - 3)(ro+1) s [ro -i,i]M + . . . (21) 

where the sum goes over all hook diagrams of the size m. However, on the Jones topological locus, 

Ak _ A—k , 2 

P>=Pl = ^=iT ^ (22) 
only ST , = % \ - and ST , , n = ST 01 = - are non- vanishing so that 

[n] {<?} [n-2] {q} b 



„-(m+l)(m+2) . s 
H lm,m+1]( A = q 2^ = A -(m 2 -l) H [m,m+l]| p * } = « _ ^4m+4 _ q 2m+2 _ q 2m + g 2^ = 



-(m+l)(m+2) 

= = J (m) (23) 

q-q- 1 

The four-term polynomial J^ m \q) satisfies the 4-term recurrent relation 

j(m) _ q i j(m-l) _|_ q -2m+8 ^ j(m-2) _ j(rn-3)^j ^4) 

For generalizations to Jones unreduced superpolynomials for the same series of knots see [17] (where slightly 
different normalization is used and q — >• \jq\ hereafter, such a difference with expressions for the trefoil in the 
literature arises, because we use the right-hand trefoil described by braid (1, 1|1, 1), while often the left-hand 
trefoil (—1,-1 1 — 1,-1) is used instead). 
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In fact, there is nothing special about the [to, m + 1] series, at least, for the HOMFLY polynomials. For 
[to, to + r] one similarly gets 



/ X -(m+l)(m+l+r) 
jj{m,m+r\t^ _ 2| \ _ -2(m-l)(ro+r) / (ro-l)(m+r) g* _ (m-3)(m+r) g* \ _ J_ j{m\r) ^5) 

n v y \hi h \h [m] y [m— 2] y 1 v ; 



-(m+l)(m+l+r) 

where the 4-term (differently normalized) Jones polynomial 

j(m\r) _ q 4m+2r+2 _ q 2rn+2r _ q 2m + q 2 ^6) 

satisfies 

rj-{m\r) q 4r rj-(rn—r\r) q —2m+2r+6 ^ >j{ra — 2\r) rj- (m— 2 — r|r)^ (27) 

Note that this relation at some values of to and r involves both the Jones polynomials of knots and of links. For 
instance, at to = 7 and r = 3, T[7, 10], T[4, 7], T[5, 8] and T[2, 5] are involved, all of them being knots, while at 
m = 8 and r = 3, of involved T[8, 11], T[5,8], T[6,9] and T[3,6] only the first two are knots, while the second 
two are links. 

3.3 Colored HOMFLY by cabling 

Another application of (7) and (17) is to construct colored HOMFLY polynomials from the ordinary ones. To 
this end, as a first step, one has to use the cabling procedure, i.e. the relation 

Hft® m — Hft (28) 

where JC m denotes m-cabling of the knot K. (see, e.g., [18, eq.(1.5b)]). At the same time, the l.h.s. of this 
formula can be presented as a sum over the irreducible representations, 

H^ m =^H§ (29) 

s 

Therefore, in order to calculate Hg, one has to calculate JC m in the fundamental representation and then project 
out the result onto an irreducible representation Q. This is done using proper projectors. 
For instance, relation (7) implies that 

p 1 + q7Z a p <f-qR 

[2] = TT7" [11] = TT?~ ( } 

are projectors correspondingly onto representations [2] and [11], which emerge in the decomposition of the 
product [1] x [1]. Similarly, (17) implies that the projectors onto irreducible components of the triple product 
[1] x [1] x [1] = [3] +2[21] + [111] are 

_ (i-g 2 n 3 )(nl + n 3 + i) _ e (n- q 2 )(nl + n 3 + 1) 

" (l-^Xl+^+g 4 ) ' (l_ g 4 )(l+g2+ g4 ) W 

and 

p± _ , {q 2 U3 1){n3 _ q 2 ){Uz _ e ±2W3) 

[21] — ^ l + q* + <? {62) 

projects on two representations [21] that correspond to the two eigenvalues e ±27 ™/ 3 of 1Z 3 . 

4 Differentials 

Another type of relations comes from the trivial observation that 

H*(A\q) H%(q N \q) = (A q N )G^(A\q), 
H*(A\q) H%(-(l/q) N \q) = (A- {-l/q) N )G*(A\q) (33) 
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with some polynomials G^(A\q), G^(A\q). This naturally leads to the formalism of " differentials" , widely used 
in application to superpolynomials. It is effective, when an additional information is used about the HOMFLY 
polynomials for a particular N, for example, that Sq = Sq{pI} in (4) arc non-vanishing only for l(Q) < N (in 
particular, for TV — 1 a single term with Q — [m\R\] survives in the character expansion of H^(q\q)). 

These relations have a nice pictorial representation [15]. Let us plot the HOMFLY polynomials H^(A\q) as 
a set of points corresponding to the monomials A k q l on the (k, l)-plane with positive and negative signs denoted 
by the white and black circles respectively, and with multiplicities (if necessary). Then, all the terms that could 
cancel under the specialization A — q N , lie on the lines Nk = — I + const, and they cancel if the algebraic sum 
of the terms on each line vanish. Similarly, all the terms that cancel under the specialization A = —l/q N lie on 
the lines Nk = I + const, and the algebraic sum of the terms on each such line should also vanish. An important 
point is that this line structure is preserved by a /^-deformation, and this provides an important tool to deal 
with the superpolynomials [15]. 

As the simplest example, the fundamental HOMFLY polynomials for all the 2-strand knots pictorially look 
basically the same: since 



H& 2n+1] (A\q) q 2n + 2 



Q 



S* 



-A — 



(34) 



they are sets of n + 1 and n points on the horizontal lines k — +1 and k = — 1 respectively, and one has, say, 



for the [2, 9] case, for (q s 



q 4 + l 



+ q~ 8 )A - (q 6 + q 2 + q- 2 + q~ 6 )A- 




A = q (N = 1) 




A = -\ (N = l) 



A k t 



o o 



A = q 3 (Y = 3) 





A = q 7 (N = 7) 



The circles not connected by oblique lines survive to describe the specialized polynomials H n (A = q N \q). 
even N and for odd N > 9 no reduction takes place: all the nine items are present. 



For 
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Similarly, for the figure eight knot 4i with ifC 1 -- 1 ! 1 -- 1 ) = 1 + {Aq}{A/q} = A 2 - q 2 + 1 - q~ 2 + A~ 2 one has 
A k i 




A = q 
A = -l/q 



(N=l) 



and there are no non-trivial reductions, except for A = q and A = — 1/q (N = 1). 

One can also consider "unreduced" HOMFLY "polynomials", not divided by Sr. If A is considered as an 
independent variable, they are actually series, not polynomials in q and their pictorial representation contains 
(semi)inhnitely many vertices, of which finitely many survive after specializations A = q N with arbitrary 
N = 1, 2, 3, . . . For example, already for the trefoil 3i one has 



unreduced 



H^(A\q) = -q^ T A ~ 




2 +q- 2 )A-A- 1 ) 

j 

V 

reduced 



(35) 



and 





A = -1/q 



Here we manifestly indicated non-zero sums £ along the lines, which, hence, survive after the specialization. 

Deformation to superpolynomials dictated by these pictures by the rules of [15] is actually different for the 
reduced and unreduced HOMFLY polynomials, what provides two different families of superpolynomials: the 
reduced and unreduced ones. The former are much simpler and have more profound structures (in particular, 
possess character decompositions and can be extended to arbitrary time variables [10, 6]), however, the latter 
seem to be related to somewhat simpler versions of Khovanov-Rozhansky homological descriptions. Of course, 
the two versions are not independent, for the same example of the trefoil in the fundamental representation one 
has 



H 3 a 1 (q) = { 



reduced : P* 1 (a, q, t) = 1 + a 2 tq 2 + a 2 q 2 t 3 

unreduced : P Sl (a, q, t) = - — — - a 2 q 2 t 3 + fat + -) - l) 

q q^ 1 V a/ \q^ / 



(36) 



and 



^□ 3l (a,q,t) = 



a — a 



— q 2 a 2 t 3 + 



a + a H 1 



1 



q 2 a~ 2 



T (P 3l (a,q,t)-q 2 a 2 t 3 



(37) 



q — q * q — q 1 1 + q 2 a 2 t 

One can see that the unreduced superpolynomial P 31 (a, q, t), in contrast with the reduced one, contains 
non-positive coefficients and only after specialization a = q^ becomes a polynomial with positive coefficients. 
For instance, at a = q 2 (Jones) 



P 3l (a,q,t) =q 3 t + i+q 7 t 3 



q 5 t 3 



(38) 
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which, indeed, coincides with the unreduced Khovanov homology [19]. 

A relation similar to (37) persists for other knots. For instance, for any "thin" knot [15] 

*<-*'>-h£ (r +^t^&>^> -©* ) <»> 

where S/e is the signature of the knot (for the torus knot T[m, n], e.g., S/c — (m — l)(n — 1)). To make formula 
simpler, we used here slightly different normalization of superpolynomials. 

5 Pliicker relations (KP integrability) 

The third important kind of equations is related to integrability. For example, for the generating function of 
the extended HOMFLY polynomials 

{Y,B^RU\Pk}SR{pk} for knots 

(40) 
{#k}IlU SrM } } for links 

to be a r-function of the KP hierarchy in the pk //e-variables [20], the coefficients £q = Yl< rCrq^S R{q\pk} 
should satisfy the infinite set of quadratic Pliicker relations: 

£[22]Z2[0] - £[21]£[1] + C[2]£[ll] = 0, 
£[32] £[0] - £[31] £[1] + £[3]£[11] = 0, 
£[221]£[0] - £[211]£[1] + £[2]£[111] = 0, 

(41) 

This indeed happens for one special class [6]: the torus knots/links with all ay = 1 (1 < i < m — 1, 1 < j < n) 
but not in general. 

Another possibility is that "H^ themselves satisfy the Pliicker relations as functions of R, i.e. that Z® is a r- 
function w.r.t. times pk/k. Then, in particular, the Ooguri-Vafa partition function H^{q\A}SR{pk} would 
be a KP r-function in these time variables. Again, this is not true in general, but an example is known when 
this is true (besides the trivial example of the unknot): in the limit q — > 1, when the HOMFLY polynomials 
turn into the "special" polynomials 

^=«sHm (42) 

which have a very simple dependence on R: 

6%(A) = (e*(A)) m (43) 
The full HOMFLY polynomial is restored by action of a combination of the cut-and-join operators [21] 

J2H^(A)S R {Pk} =ex P [^{q-q-'r-^^WQiPk}] cxp (44) 

r \q K(a)) J \ k feV ' ) 

with 

a^{A) = &n(A) (45) 

and resembles very much the Hurwitz partition function [22, 23], which is also in general not a KP/Toda r- 
function in the variables {p,p}. In the Hurwitz case the way is known to cure this problem [24] (by switching to 
times associated with a special basis in the space of cut-and-join operators), perhaps, this analogy can show a 
way to proceed in the case of HOMFLY polynomials as well. The integer coefficients of the special polynomials 
o-q(A) are related to the Ooguri-Vafa numbers [25, 8], but in a somewhat complicated way for R ^ □. 
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6 .4-polynomials 

The fourth, currently most interesting type of relations arc difference equations in R [26], so far found empirically 
mostly for symmetric representations R — [p] [27, 28, 29] or for the small groups with N — 2, 3 [30, 31]. 

The usual strategy to derive these equations [30] (at least at TV — 2, when R = [r — 1] is given by the single 
positive integer) is to 

• convert the Jones polynomial (with A = q 2 ) into a form of a proper g-hypergeometric series [32, 33] 

M«) = e n !5' 9 2? b+(r,k) g c(r ' k)Yk ( 46 ) 

where b^ 1 (r, k) and C (r, k) are linear forms and a quadratic form accordingly, which, as well as B-f and the 
r-vector Y, are just parameters which describe the concrete knot and the g-Pochhammer symbol 

r-1 

{B,q) r = '[[{l-Bq i ) (47) 

• write down a difference equation, which always exists for the proper g-hypergeometric series. 
For example, in the simplest case of the trefoil (3i knot) [35] 

oo 

J^(q) = q 2 - 2r J2^ 2kr (^ 2r ^ 2 ) . ( 48 ) 



k=0 



the equation is [30, 34]: 

J^(q)-U r (q)4l 1 (q) = V r (q), 

1 - n 2r ~ 2 1 — n 4r - 2 

UM = -< l i - &A - Y ^-, Vr( q )= q ^ 1 - T ^ r (49) 

For arbitrary torus knot [m, n] (trefoil is either [2, 3] or [3,2]), the equation is of the second order, however, for 
the whole scries T[2, 2s + 1] [36] 

jj^Mfci+i-ar) I ' Jk ^ I (so) 

( q2 > q2 ) k \ q2 > q2 ) ki+1 - k J 

it reduces to the first order and looks similar to (49) with [36] 

U r (q) = - g 2(.+l)-2(2- + l)r y Aq) = ^. + l)(l-r) 

For non-torus knots the formulas get longer, and it is convenient to write them in terms of the operator M : 
Mf(q, r) = q- 2r .f(q, r) and the shift operator L: Lf(q, r) = f(q, r + 1), LM = q 2 ML. Then, for knot 4i 

OO 

j> (5) = £m) v fe(fe+1) ( q 2 - 2r , e) k (V +2r , e) k (52) 

fc=0 

the equation is [30, 34]: 

A 4l (L,M)J^(q) =B(M), (53) 

A = q 4 M 2 (l - M)(l - q 6 M 2 ) - (q 2 M + 1)(1 - q 2 M - q 2 M 2 - q 6 M 2 - q 6 M 3 + g 8 M 4 )(l - q 2 M) 2 L+ 

+q 4 M 2 (l - q 2 M 2 )(l - q 4 M)L 2 , B = q 2 M(l - g 6 M 2 )(l - q 2 M 2 )(l + q 2 M) 

In the limit of q — > 1, A(L, M) can be considered as an ordinary c-number function A(l,m) of L — > Z and 
M — >• m and coincides (this is called the AJ-conjecture [26]) with the ,4-polynomial [37] (this is why A K (L, M) 
is called the non-commutative ,4-polynomial), which defines a "spectral curve", associated with the knot /C: 

Y, K : A K (l,m)=0 (54) 



9 



which is naturally equipped with the Seiberg-Witten differential 

dS K = logZ d log to (55) 

According to the standard dictionary of the Seiberg-Witten theory [38, 39, 40], this form of the SW differential 
is a clear sign of the relation to 5d YM theories. 

7 The large- 1 R | limit 

Given a difference equation (53), a natural question to ask is about various kinds of the loop (quasiclassical) 
expansions for its solutions. Loop expansion arises when one introduces the 't Hooft coupling constant u by 
putting q = e h = e 9 2 = e u ^ and \R\ = \ = £ . 

7.1 Special polynomial expansion 

To begin with one can look at the anzatz, provided by the special polynomial expansion (44). It would imply 
that the solution behaves like 

Wg~o5(i4)l fl l =exp(|iZ|.loga£(A)), (56) 

i.e. grows exponentially with \R\, in accordance with the volume conjecture [41]. However, things are not so 
simple and strongly depends on the range of values of u. 

First of all, corrections to this formula could also contribute to the exponential growth. Indeed, the next 
correction is [21] 

where ~ = ^- at large \R\ = r and e = q — q^ 1 . Since the latter is equal to 2sinh?i = 2sinhg 2 ~ r^r, 



ua*~ (A) 

the exponent in this formula is proportional to ckr <~ \R\ at large \R\, so that one should add to 

log o~n(A) in (56). The next correction will be proportional to —. — [3] , a and so on: there is an infinite series 

in powers of u of contributions to the leading (<~ |i?|) "volume" behavior of the exponent. 

Second, for N fixed, i.e. for A — > 1 (for example, for the Jones polynomial with A = q 2 ), this volume 
contribution is in fact vanishing term by term, because a^(A = 1) = 1 and all other (Tq(A = 1) = [21]. In 
fact this is in accordance with the well known fact that for small enough u there is no exponential growth in 
J r (q = e"/ r ). Instead (this formula was realized for the figure eight in [42] and for generic knots in [43] basing 
on the Melvin-Morton-Rozansky conjecture [44]) 

*<■ - ^ - - (§ ® " H - + r { ^Z,? m 

with some polynomials Wk- In particular, 

fo(u) = - log (Alexander^ = e")j (59) 

and this is in a nice accordance with the special polynomial expansion [21]. 

Only when e u exceeds the smallest root of the Alexander polynomial, another solution appears, with /_i(u) 7^ 

0. 

At the same time, if one chooses u = 2iri, there is the exponential behaviour (56) [45] and the resulting 
coefficient in front of \R\, for the Jones polynomial A — q 2 , is equal to the hyperbolic volume of the knot [41] 
(volume conjecture). 
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7.2 Loop expansion 



It is described by the standard loop expansion: 

J r (q = e<?) = cxp ^ .9 2 "A(«)j (60) 
where v = u — 2ni, i.e. v = is the volume conjecture point. The shift operations are 



J r±1 {q = e g2 ) = cxp 9 2k fk{v ± g 2 )\ = exp 



,A;>-1 



£ {-yg 2k+2] dlh(v) (61) 



\ 3 = 

\fc>-l 



Therefore, the difference equation (53), if expanded in powers of g 2 , turns into an infinite system of recurrent 
relations for infinite set of functions fk{v)- Remarkably this system can be related [46] to a similar system of 
equations for multidensities p p ,k{vi, ■ ■ ■ , vu) 1 which constitutes the essence of the AMM/EO topological recursion 

[47]. More exactly, the knot polynomial is associated with ^ log det(M— u)^, where the brackets denote averaging 
over the "matrix" M in the would-be matrix model underlying the recursion. 

Recursion itself means that, given /_i(w), usually described in terms of the spectral Riemann surface (54) 
with a Sciberg-Witten differential (55), one can recursively reconstruct all the functions fk{v). 



8 M/M AGT relation 

8.1 ^4-polynomials, state integral model and Baxter equations 

As soon as /_i ^ appears, one can neglect the r.h.s. in equation (53) and substitute it by the homogeneous 
.4-polynomial equation. This one has a simpler solution, in the form of integrals of the double-sine ratios instead 
of sums of the q-Pochhammer ratios [48] . This solution is known as partition function of the state integral model 
[49]. 

According to [34] one can interpret this homogeneous difference equation as a Baxter equation for some 
relativistic integrable system (such systems are associated through Seiberg-Witten theory [50, 38, 39, 40] to 5d 
SYM theories, what explains the possible name for such a relation). 

In fact, this should not come as a surprise due to the well-known connection between 2d CFT and Chern- 
Simons theory. As soon as the former one is related via the AGT correspondence [51, 52, 53] with the 4d SW 
systems, one would naturally expect a similar relation in the dimensions increased by unit: the one between 
Chern-Simons theory and 5d SW theory 2 . On the other hand, the Baxter equation is an evident object to look 
at, since it encodes changing traces of the group element with representation (fusion), and there is a counterpart 
of this on the knot theory side: the difference ,4-polynomial equation describes exactly changing the traces (of 
the braid group element) with representation. The technical similarity of these two things is however appealing. 

To be more concrete, one can notice that eq.(53) has the form of the Baxter equation for the relativistic 
Toda chain, which, in accordance with [39], describes the 5d SW system. Instead, one could definitely start 
from the classical limit of the Baxter equation (relativistic Toda spectral curve) noticing that it coincides with 
the "classical" ,4-polynomial, and then apply the topological recursion with the 5d SW differential dSsw (55). 
This was done in [46] and the result reproduced the expected asymptotics of the Jones polynomial (see also [55] 
for a more generic context). 

This is depicted at the scheme: 

AH = — > bd Baxter(Q) = 

dS = d\og^! 



t H = Q I h = (62) 



H AMM/EO topological recursion E : A = 

< — dSsw 



2 Within a different context, a 5d/3d duality emerged in [54]. 
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This scheme implies that, at the quantum level, there is a relativistic integrable system/5d SW system 
described by the Baxter equation for the Baxter operator Q(fi). The Fourier transform of it, ^(x) satisfies a 
relativistic Schrodinger equation, and one can construct the 5d Nekrasov functions (in the Nekrasov-Shatashvili 
limit [56]) as the Bohr-Sommerfeld integrals of this equation [57], which are A- and B-periods of the differential 
dS = dlog^(x). On the other hand, the same Q(/z) is the partition function of the state integral model, since 
this latter satisfies the same difference equation. 

At the classical level, the system is described by the Riemann surface (the spectral curve of the classical 

integrable system) of the 5d SW system with the differential given by the classical limit of dS: dS ^> dSsw- 
This Riemann surface coincides with the classical *4-polynomial (which describes the SL(2) representation of the 
knot complement to S* 3 , as viewed from the boundary [37]) and dSsw (55) allows one to apply the topological 
recursion and to restore the partition function of the state integral model 

8.2 Universal group element 

Let us now briefly describe the origin of the Baxter equation as the fusion relation of the universal group element. 
First of all, one has to construct the group element given over the non-commutative ring, the algebra of functions 
on the quantum group (see a review in [58]). That is, one constructs such an element g £ U q (Q) ® A(G) of the 
tensor product of the Universal Enveloping Algebra (UEA) U q {Q) and its dual algebra of functions A{G) that 

&u{g) = g®ug& A{G)®u q (g)®u q (g) (63) 

To construct this element [59, 58], we fix some basis in U q {G). There exists a non-degenerated pairing 
between U q (Q) and A(G), which we denote < ... >. We also fix the basis in A(G) orthogonal to 

w.r.t. this pairing. Then, the sum 

T = J2 x(a) ® Tia) ^ A (G)®U q (g) (64) 

a 

is exactly the group element we are looking for. It is called the universal T-matrix (as it is intertwined by the 
universal 1Z- matrix) or the universal group element. 

In order to prove that (64) satisfies formula (63) one should note that the matrices M" 13 and giving 
respectively the multiplication and co- multiplication in U q (Q) 

T (a) , T ((3) = M a0 T ( 7 ) } A(T< a ') = D^T^ ® (65) 
give rise to, inversely, co-multiplication and multiplication in the dual algebra ^4(G): 

Df^ = (A(T( a l),J:M0JfW) = (t("U M ^ (7) ) (66) 
Mf = (t^T^\X^ = (tW 8rW,A(JfW)) 

Then, 

Au(T) = x(a) ® &u{T {a) ) = ^ D h x[a) ® T(f>) ® T<7) = Jf w JT (7) ® T {f>) O = T ®u T (67) 

This is the first defining property of the universal T-operator, which coincides with the classical one. The 
second property that allows one to consider T as an element of the " true" group is the group multiplication law 
g ■ g' = g" given by the map: 

g ■ g' = T <g) A T e A(G) ® A(G) <g> U q (G) — > g" e A(G) ® U q {Q) (68) 

This map is canonically given by the co-multiplication and is again the universal T-operator: 

T ® A T = x(a) ® x(l3) ® T {a) T i(1) = M l,a x(a) ® xi0) ® T<7) = J2 A ( x(a) ) ® T (a) (69) 

i.e. 

5 = T(A,T), ,g' = T(A',T), 5 " = T(A",T) 
lE{l(«'®/}e4(G)®J, A' = {J® A^} e J<8> A(G) I"E{A(l(°))}e4(G)®i(G) (70) 
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8.3 Baxter equation 

The most essential property of the group elements is that they are intertwined by the universal i?-matrix [60] , 
i.e. 

K T® V T' = T'o c/ T Tl (71) 
This means that their traces taken at any representations of the UEA, Ty commute: 

[T V ,T V ,]=0 (72) 

The equation that connects these traces in different representation is called fusion and its specific limit is 
the Baxter equation. As a simple instance, let us consider the case of quantum affine algebra SL(2) and the 
evaluation representation of this latter Vj(X) labeled by the spin j and evaluation parameter A (see, e.g., [61] 
for a review) 3 . The product of two such representations (j, A) and (k, fj,) is an irreducible representation unless 

- = qi+ k -P+ 1 (73) 
A* 

where p = 1, min(2j, 2k). In this latter case, the product is partly reducible: 

o — ► v^ p/2 ( q p/ 2 x) ® T4- p/2 (<r p/ V) — ► Vj(\) ® v k {fj) — ► 

— ► V} +fe - p/2+ i/ 2 (<z p/2 - fc - 1/2 A) ® ^/ 2 -i/2(5 fe+1/2 - p/2 M) — ► (74) 
In terms of traces of the group element, this (fusion) relation can be written in the form 

Tj(A)Tfc(/x) = T j - p/2 ( g f/ 2 A)T fe _ p/2 (g-f/ 2 /i ) + T, +fe _ p/2+1/2 ( <Z f/ 2 - fe - 1 /2A)T p/2 _ 1/2 (g fe + 1 /2-p/2 M) (75) 
Choosing fc = 1/2, p = 1 and A = qi +1 / 2 (j,, one immediately obtains 

T^+^T^) = T j _ 1/2 (^+ 1 /U ) + T j+1/2 (g>) (76) 
It turns out [62] that (upon the proper normalization of Tj(A)) there is a finite limit 

lim ■W+V^) = Q(/x) (77) 

and (76) turns into the standard Baxter equation 

Ti/ 2 (/i)0(/i) = 0(«A*) + (78) 

Here T 1 / 2 (/x) is a Laurent polynomial of degree 1 both in ji and fi . One can generate Laurent polynomials 
of higher degrees considering in (76) traces of products of the group elements, (68), the product of n group 
elements giving rise to degree n. Thus, generically one would get instead of (78) 

Pj ft ^ 1 ]Q(/i)=Q(p) + Q(fV) (79) 

which coincides with the Baxter equation for the n-particle relativistic Toda chain, which corresponds at the 
SW level to the 5d pure gauge theory. The transition to four dimensions and the standard Toda system is done 
by replacing the affine algebra with the Yangian. This leads to the same equation with P n being a degree n 
polynomial of \i. 

8.4 Knots, Baxter equations and 5d SW theory 

One can also choose p = 2 and k = 3/2. Then, 



Q0i)T 3/2 (/i) = t 1/2 (<tV)Q(<?>) + T 1/2 ( w )Q(<rV) 



(80) 



This equation is another version of the Baxter equation, it gives rise to the same spectral curve in the classical 
limit. 



3 Remind that these representations can be treated as either finite dimensional but reducible, or infinite dimensional irreducible 
representations [61]. We prefer here the latter point of view. 
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Let us consider, for instance, the first non-trivial example of the figure eight knot 4i which ,4-polynomial 
equation is the homogeneous part of (53), i.e. 

A Al {L,M)Z Al {M) = (81) 

We use the basis where M acts on 2 diagonally so that (81) reduces to 

(i - «V)«(<rV) + (i - sVWri = ftke-'lOM , 

^-■i^v+^H^+^M^+f) (82 » 

where Q{p) = (1 — /i)3 4i (m) with identification M = \i. (82) is nothing but (80) with 



Ti/ 2 (/i) = ? 2 M - -i- , T 3/2 (/i) = ( q 2 fi - -L) P 2 \n, [T 



(83) 



The coefficients of T 1 / 2 (m) and T 3 / 2 (m) can be related with concrete values of the energies (integrals of motion) 
in the periodic 2-particle relativistic Toda chain, and the corresponding SW curve (5d pure gauge theory) is 

i+ H* 2 v)-H)- 2 < M > 

In contrast with the generic SW curve, where there is an arbitrary degree two Laurent polynomial at the r.h.s. 
of (84), a peculiar polynomial with all coefficients fixed corresponds to the figure eight knot. Similarly, (83) 
specifies very concrete polynomials, while in the generic integrable system these are arbitrary ones (of fixed 
degree). 

Note that choosing k = 1, p = 1 and A = qi +1 /x (or p = 2 and A = q^fj,), one obtains from (75) a third order 
difference equation (or its counterpart with q — > l/q) 

Titf^MM) = <2(<? 2 M) + T 1/2 (g M )Q(<rV) (85) 

The procedure can be continued further to generate higher order difference equations which can be compared 
with the ,4-polynomial equations for other, more complicated knots. In accordance with [30], such an equation 
exists for any knot. The coefficients of the polynomial are fixed by choosing the representation of 

the algebra of functions on the quantum group. Indeed, due to (72) we know that all the coefficients of this 
polynomial are commuting. They are nothing but the Casimir functions that determine the representation of 
A(G)._ 

In integrable systems these coefficients are just integrals of motion, which can have arbitrary values. On the 
contrary, in knot theory their values are fixed to be very concrete numbers, i.e. the representations of A(G) in 
knot theory are fixed. It still remains unclear what is the condition that fixes them in group theory. 

Now let us stress out that all these Baxter equations or quantum ,4-polynomials correspond in the classical 
limit to 5d SW theories, generically of the quiver type [63, 64]. On the integrable side, these theories are 
described by various degenerations of the sl(p) XXZ spin chain [40] (or corresponding Gaudin type systems, 
see [65]). Their quantization is immediate, either as integrable systems which goes exactly through the Baxter 
equation [57], or as ,4-polynomials [46, 55]. 

Even more general Baxter equations are described by the generic sl{p) XXZ spin chain, in the classical limit 
giving rise to the 5d SW curves that are plane complex curves given by sets of arbitrary Laurent polynomials 

PnA^V- 1 ] [40] 

P 

^VPn^/i- 1 ] =0 (86) 
i=l 

However, obtaining these Baxter equations directly from the representation theory as in the previous subsection 
is not that immediate in this case, and the whole construction becomes less elegant and much more involved 

[66]. 

Note that one can similarly consider the Baxter equation describing SL(N) . In terms of knot theory, it 
should describe the difference equations for the HOMFLY polynomials with specialization A — q N . 



14 



8.5 State integral model and open relativistic Toda chain 

Let us note that solutions to the Baxter equations are usually quite involved, and are unknown in a generic case. 
On the contrary, the partition function of the state integral model, which solves the *4-polynomial (i.e. Baxter) 
equation is often known. The reason is that with peculiar values of the integral of motions, i.e. with the specific 
form of traces T, corresponding to knots, the equation is explicitly solved. The answer is presented as a multiple 
integral of product of the quantum dilogarithms [67]. This answer is in fact associated with the wave function 
of the open relativistic Toda chain (which, within the SW context, corresponds to the perturbative limit of the 
gauge theory). More exactly, this wave function solves the Schrodinger equation given by the Hamiltonian of 
the system dual [68, 69, 70] to the relativistic Toda chain. In other words, this is the (difference) equation that 
describes the dependence of the wave function on the energies. 

It happens, since the Baxter equation, which is a difference equation for the separated variables in the 
periodic chain, is simultaneously (at peculiar values of energies) a Schodinger equation for the (dual) open chain 
(this is possible definitely only due to the fact that chain is relativistic and, hence, both the Schrodinger and 
the Baxter equations are difference ones). 

Let us see how it works in already discussed simplest non-trivial example of the figure eight knot. Remind 
that the Hamiltonian of the fc-particle relativistic Toda chain is given by [71] 

k 

h{ Xi ; Pi } = [ e " m + q- 1 g 2 ^e^- x ^ ]+1 -^ ul ^ + ^ p ^ (87) 
i=l 

with the boundary condition Xk+i — X\ in the periodic case and Xk+i — oo in the open case. Here q = e™ 1 ^ 2 
is the relativistic parameter (q — > 1 as the speed of light goes to infinity), e = 0, ±1 parameterizes different 
relativistic Toda systems and g is a coupling constant, which can be removed by redefinitions of coordinates Xi 
only in the open case. In the 2-particle open case, the Schrodinger equation for e = can be reduced (upon 
neglecting the [/(l)-factor) to the equation 

V 7 (<rV) + V> 7 M + 9 2U1 ^(P) = (e" 7 + e^ 7 ) iM/x) (88) 

where 7 parameterizes the energy and we fixed u 2 = 1, oj\ = h/iri and denoted (1 = exp(27rx), x = X\ — x 2 . A 
solution to this eigenvalue problem is of the form [67] 



^> 7 (/x) oc / dt 



c 







log g 


7TZ ' J 


s 


[it + q-^siogg 


A,l) 

7TI 5 J 



7T 2 T (2t + x) 

e R (89) 



where S(x\h/wi, 1) [67] is related with the quantum dilogarithm [49] by the formula 



. . „ / z iri + h n \ 



and the contour C is defined in [67, (2.26)]. Now, if one puts ng 2uJl — 1, integral (89) reduces to (^ = 

V 7 (<T 2W1 ) oc / dp $ft(p + frr + ft) ^ {p+7n/2) _ vl/2 

which is exactly the partition function of the state model for the figure eight knot as a function of u = ttj/2 
[49] 



a 4i (e M )cxv> 2tlA (. 9 - 2wi ) 



(92) 



Since the dependence on 7 is regulated by the dual Toda Hamiltonian, we come to the claim above. 

In the case of more complicated knots one has instead of (89) a multiple integral of the product of the 
quantum dilogarithms [49]. This integral could be associated with multi-particle relativistic open Toda wave 
functions [67] and their further generalizations. 
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9 Conclusion 



To conclude, in development of the program outlined in [34, 6], we presented a brief review of the new and 
very interesting research field: the study of interrelations between various HOMFLY polynomials, probably 
reflecting the powerful underlying integrable structure, which still remains to be fully revealed. The relations 
involve different knots/links and different representations, they arc rather diverse and form a somewhat strange 
intertwined set, and partly unexpected parallels are seen with the Hurwitz theory (i.e. with that of symmetric 
group characters), quantum integrable systems and AGT relations. Even more interesting should be general- 
izations from the HOMFLY to swperpolynomials [14], however, what structures survive this /3-deformation [72] 
remains an open question. 

A way to systematical study of these problems is now open by the powerful HOMFLY calculus, developed 
in [10, 6, 7, 8, 9], exploiting the simple and universal structure of SU q i?-matrices and Racah coefficients, and 
allowing one to obtain the parameter-dependent expressions for knot polynomials. One now expects a fast and 
fruitful development of this old and intriguing field. 
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